Fluid dynamics and jamming in a dilatant fluid 
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We present a phenomenological fluid dynamics model for a dilatant fluid, i.e. a severe shear 
thickening fluid, by introducing a state variable. The Navier-Stokes equation is coupled with the 
state variable field, which evolves in response to the local shear stress as the fluid is sheared. The 
viscosity is assumed to depend upon the state variable and to diverge at a certain value due to 
jamming. We demonstrate that the coupling of the fluid dynamics with the shear thickening leads 
to an oscillatory instability in the shear flow. The model also shows a peculiar response of the fluid 
to a strong external impact. 

PACS numbers: 83.80.Hj,83.60.Rs, 83.10.Ff,83.60.Wc 



Dense mixture of starch and water is an ideal mate- 
rial to demonstrate the shear thickening property of the 
non-Newtonian fluid. It may behave as liquid but is im- 
mediately solidified upon sudden application of stress, 
thus one can even run over the pool filled with the fluid. 
It is amazing to see that it develops protrusions sticking 
out of the surface when it is subject to strong vertical 
vibration [J HI) and also intriguing that the fluid vibrates 
spontaneously when one simply pours it out of a con- 
tainer. A source of these unusual behaviors is the severe 
shear thickening; the viscosity increases almost discontin- 
uously by orders of magnitude at a certain critical shear 
rateQ, which makes the fluid behave like a solid upon 
abrupt deformation. Such severe shear thickening is of- 
ten found in dense colloid suspensions and granule-fluid 
mixtures 043 j and they have been sometimes called "di- 
latant fluid" due to the apparent analogy to Reynolds 
dilatancy of granular media[8|. 

Despite its suggestive name, physicists have not 
reached the microscopic understanding of this shear 
thickening property. It was originally proposed that the 
shear thickening is due to the order-disorder transition 
of dispersed particles 0, H [l(|; the viscosity increases 
when the fluid flow in high shear regime destroys the 
layered structure that has appeared in lower shear rate 
regime. Although there seem to be some situations where 
this mechanism was believed to cause the shear thicken- 
ing, there are other cases where the layered structure is 
not observed prior to the shear thickening pdj or no sig- 
nificant changes in the particle ordering are found upon 
the discontinuous thickening [l^, EH- The cluster for- 
mation due to hydrodynamic interaction [l4l - fl7| and the 
jamming [1, [f| [l8l - [22T | have been proposed as plausible 
mechanisms. 

There are several peculiar features in the shear thick- 
ening shown by a dilatant fluid: (i) the thickening is so 
severe and instantaneous that it might be used even to 
make a body armor to stop a bullet |23|, (ii) the relax- 
ation after removal of the external stress occurs within 
a few seconds, that is not very slow but not as instan- 



taneous as in the thickening process, (iii) the thickened 
state is almost like a solid and does not allow much elastic 
deformation unlike a visco-elastic material, (iv) the vis- 
cosity shows hysteresis upon changing the shear rateQ, 
(v) noisy fluctuations have been observed in response to 
an external shear stress in the thickening regime fA |22|. 

In this report, we construct a phenomenological model 
in the macroscopic level and examine the fluid dynamical 
behavior of the dilatant fluid. We introduce a state vari- 
able that describes phcnomcnologically an internal state 
of the fluid, and couple its dynamics with the fluid dy- 
namics. We find that the following two aspects of the 
model are important: (1) the fluid state changes in re- 
sponse to the shear stress, (2) the state variable changes 
with the rate proportional to the shear rate. We examine 
the model behavior in simple configurations and demon- 
strate that the model is capable of describing the char- 
acteristic features of hydrodynamic behavior and shows 
the oscillatory instability. 

a. Model: The model is based upon the incompress- 
ible Navier-Stokes equation with the viscosity ?y(</>) that 
depends upon the state variable <fi. The scaler field <f>, 
which takes a value in [0, 1], represents the local state of 
the medium. We assume two limiting states: the low vis- 
cosity state {4> ~ 0) and the high viscosity state (cf> ~ 1). 
At 4> = 1, the system is supposed to be jammed and the 
viscosity diverges. 

The state variable cf> relaxes to a steady value 0* deter- 
mined by the local shear stress S; 4>* is supposed to be 
a continuous and monotonically increasing from to (j>M 
as a function of S with a characteristic stress So- The 
limiting value <Pm represents the state of the medium in 
the high stress limit and should depend upon the medium 
properties such as the packing fraction of the dispersed 
granules. In the following, we employ the forms 
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with a dimensionless constant A. We assume the Vogel- 
Fulcher type strong divergence in rj(<j)) in oder to repre- 
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FIG. 1: Simple shear flow configurations with the coordinate 
system. 





sent the severe thickening, but a detailed functional form 
is rather arbitrary. 

For the simple shear flow configuration as in FigQTa) 
with the flow field v given by (u(z, t), 0, 0), the model is 
described by the following set of equations, 



du(z, t) 
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d<Kz,t) 
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S(z,t), 
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~(<f>(z,t)-MS(z,t))), (3) 



with p and r being the medium density and the relaxation 
time. The shear stress S(z, t) is given by 



S(z,t) = T]((f))j, i=-^- 



du 
0? 



(4) 



where the shear rate is denoted by 7. 

Now, we suppose the relaxation of <p toward is 
driven by the shear deformation, thus the relaxation rate 
1 /t is not constant but proportional to the absolute value 
of the shear rate, i.e. 



1/t = \y\/r 



(5) 



with a dimensionless parameter r. Note that, with this 
form of relaxation, the state variable (f> does not exceed 
1 even when 4>m > 1 because 7 = at <fr = 1 due to the 
diverging viscosity. 

This particular form of relaxation is employed in order 
to represent the athermal relaxation driven by local de- 
formation. The system responds in accordance with the 
deformation rate, and does not change the state when it 
does not deform. This is natural if we consider the state 
change is caused by local configurational changes of the 
dispersed granules in the situation where Brownian mo- 
tion does not play an important role. 

b. Steady shear flow: First, we examine the model 
behavior in the steady shear flow configuration with a 
fixed external stress S c . The boundary condition is then 
given by 



So = S(z,t) 



z=±h 



(6) 



The steady state solution for Eqs.® and (0) with this 
boundary condition is obtained easily as 



(j) = <p*(S e ), 7 



So 



i*(s e ). 



(7) 



FIG. 2: The stress-shear rate relation for the viscosity given 
by Eq.(fj} with A — 1. The inset shows the same plots in the 
logarithmic scale. We employ the unit system where r/o = 
S = p = l. 



The stress-shear rate curves are plotted in Fig J2] for <jf>* 
and T](cf>) given by Eq.([T|) with A — 1 for various values of 
<j>M- In the logarithmic scale, the straight line with the 
slope 1 corresponds to the linear stress-shear rate relation 
with a constant viscosity. In the curve for cf>M = 0.8, one 
can see the two regimes: the low viscosity regime and the 
high viscosity regime. Between the two regimes, there is 
an unstable region. From this stress-shear rate curve, 
we expect there should be hysteresis upon changing the 
shear rate 7 with discontinuous jumps between the two 
branches. The jumps correspond to the discontinuous 
change of viscosity, i.e. the abrupt shear thickening. For 
the case 4>m > 1, the curves do not have an upper lin- 
ear branch because the viscosity diverges and the fluid is 
solidified. 

c. Oscillation in the shear flow: If the external 
stress S e is kept at a value in the unstable branch, where 



^<0 
dS c ' 



(8) 



i.e. the shear rate decreases as the stress increases, the 
steady flow may become unstable. The linear stabil- 
ity analysis within the solution of the form (u{z,t),Q, 0) 
shows that the mode whose wave number k in the z-axis 
satisfies 



< k 2 < p 



7* 



8S r 



(9) 



grows and the threshold mode k c oscillates at the finite 
frequency 




(10) 



Since the smallest possible wave number k min is given by 
femin = 2ir/4:h, we expect the oscillatory flow appears for 
k m in < as the system width increases if the external 
stress is in the unstable branch. 
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FIG. 3: Oscillation of the average shear rate u(h)/h in shear 
flow for h =1.3, 2, and 3 with A = 1, <f)M = 1, r — 0.1, and 
So = 1.1. The initial state is chosen as the steady solution 
© with S c = 1.0 at t = 0. 
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FIG. 4: Spatial variation of shear flow oscillation for A = 1, 
<?W = 1, r = 0.1, S c = 1.1, and h = 2. 



The oscillatory behavior of the shear flow in the unsta- 
ble regime can be seen by numerically integrating Eqs.([2]) 
and ^ with ©. In FigJ31 the average shear rate u(h)/h 
is plotted as a function of time for some values of the 
system width 2h for (f>M = 1> A = 1, and r = 0.1 with 
the external shear stress S c = 1.1 in the unstable regime. 
The initial state is chosen as the steady solution (J7J with 
S c = 1. For the case of h = 1.3, one can see the over- 
damped oscillation that converges to the steady state, 
but for the larger values of h, the system undergoes the 
oscillatory transition and the surface velocity oscillates 
with an saw-tooth like wave profile, namely, the gradual 
increase of velocity followed by a sudden drop. 

The time development of the spatial variations for the 
state variable 4> and the velocity u are plotted in Fig@]for 
h = 2. Only the positive halves (z > 0) of the symmet- 
ric solutions are shown. One can see the high shear rate 
region gradually extends towards the center as <fi relaxes, 
then the velocity drops suddenly to a very small value 
when the shear rate exceeds a certain value and <f> in- 
creases rapidly. This sudden drop of velocity is due to the 
abrupt thickening of the fluid caused by the high shear 
stress. The resulting low shear rate eventually leads to 
the low shear stress, which in turn leads to small <p with 
low viscosity, then the shear rate starts increasing again. 

d. Response to impact: The athermal relaxation ([5]) 
gives the system an interesting feature in responding to 
an external impact. We consider the simple shear flow 
but assume the fixed lower boundary at z = and the 
upper boundary at z = h is forced to move at the velocity 
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FIG. 5: The time dependence of the displacement X after 
the impact with the initial speed uq =40, 20, 10, and 5 for 
the system with <j>M — 0.8, 1, and 2. The other parameters 
are h = 2, A = 1, r = 0.1, and m = 1. 



Mo at t = (FigQJb)). Then, the motion of the upper 
wall is given by the equation, 



du(h, t) 
1 dt 



-v(<t>) £ 

oz 



(11) 



where m is the mass of the upper wall per unit length. 
Figf5] shows the displacement X of the upper wall, 



X{t) = / u{h,t')dt', 



(12) 



for the three cases, 4>m = 0.8, 1, and 2, for various initial 
speeds. The wall decelerates rapidly as the fluid thickens 
due to the imposed stress, and eventually stops. In the 
case of 4>m = 0.8, the total distance that the wall moves 
before it stops increases as the initial speed uq. On the 
other hand, in the case of cf>M = 2, the total displacement 
does not depend on it . This is because the fluid jammed 
at a certain strain as it deforms, and cannot deform fur- 
ther. 

e. Concluding remarks: There are a couple of re- 
lated models: the soft glassy rheology(SGR) model[24l| 
and the schematic mode coupling theorv(MCT) (25| . In 
the SGR model, SGR is extended to describe the shear 
thickening by introducing the stress dependent effective 
temperature, which may be compared with the inverse of 
the state variable (j> of our model. In the schematic MCT, 
the jamming transition has been examined by means of 
MCT incorporating the effects of shear schematically. 
Both are semi-empirical but deal with microscopic pro- 
cesses and produce the similar flow curves as FigJ^l On 
the other hand, the present model is phenomenological 
one only for macroscopic description to study the inter- 
play between the fluid dynamics and the shear thickening, 
which is embedded in Eqs.|T|). 

Although our model is phenomenological, the micro- 
scopic picture we may have for a dense mixture of gran- 
ules and fluid is the following; in the low stress regime 
((f) ~ 0), the particles are dispersed thus the medium 
can deform easily due to the lubrication, whereas, in the 
high stress regime (cf> ~ 1), the particles form contacts 
and force chains to support the stress, and eventually 
get jammed when the packing fraction is large enough 
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FIG. 6: Schematic illustration of configurations for low stress 
(left) and high stress(right) regime. 



(FiglHl). With this picture, it is natural to assume that 
the time scale for the state change is set by the shear rate 
in the case where the athermal deformation drives the 
configurational changes, and then the parameters Sq and 
r correspond to the stress and the strain, respectively, 
around which the neighboring granules start touching 
each other. In the case where the thermal relaxation 
takes place, the time scale may be given by 



171 
r 



1 

Tth 



(13) 



where r t h is the time scale that the thermal agitation 
dislocates the contacts. 

We have assumed that the state is determined by the 
shear stress, but it is instructive to consider what would 
happen if the state is determined by the shear rate. If 
4>„ is a function of 7 instead of Eq.JT]), then the stress- 
shear rate curve becomes monotonically increasing with- 
out an unstable branch, thus there should be no disconti- 
nuity and hysteresis, thus no instability that leads to the 
oscillation. This suggests that experimentally observed 



discontinuous transition Q and the hysteresis can be 
attributed to the "shear-stress thickening" property of 
the medium. The direct consequence of the shear-stress 
thickening may be seen in the gravitational slope flow, 
where the flux does not increases monotonically as the 
inclination angle because the fluid becomes more viscous 
under larger shear stress, thus it may flow slower as the 
slope is inclined steeper. 

A peculiar result of the present model is the oscillatory 
instability in the shear flow. Superficially, this might look 
like a stick-slip motion, but the physical origin is differ- 
ent; the stick-slip motion is caused by the slip weakening 
friction when the system is driven at a constant speed 
through an elastic device. If the present system is driven 
at a constant velocity, the system settles in one of the sta- 
ble states for a given 7. The oscillatory behavior in the 
present system appears under the drive with a constant 
shear stress, and it is a result of the coupling between the 
internal dynamics and the fluid dynamics. In this sense, 
it is also different from the oscillation in the SGR model, 
where the fluid dynamics is not considered. We found 
that the oscillatory flow also appears in the gravitational 
slope flow. 

Regarding the oscillations in experiments, one can eas- 
ily notice the vibration around 10 Hz when pouring the 
cornstarch-water mixture out of a cup. In the litera- 
ture, the pronounced fluctuations in the deformation rate 
have been reported in the shear stress controlled experi- 
ment near the critical shear rate M. l22ll and the pressure 
driven flow through microchannels[26], but not many ex- 
periments have been reported yet. 
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